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Abstract: Quantitative assessment of the impact of climate variability and human activities on runoff 
plays a pivotal role in water resource management and maintaining ecosystem integrity. This study 
considered six sub-basins in the upper reaches of the Yangtze River basin, China, to reveal the trend of 
the runoff evolution and clarify the driving factors of the changes during 1956—2020. Linear regression, 
Mann-Kendall test, and sliding +test were used to study the trend of the hydrometeorological elements, 
while cumulative distance level and ordered clustering methods were applied to identify mutation points. 
The contributions of climate change and human disturbance to runoff changes were quantitatively 
assessed using three methods, i.e., the rainfall-runoff relationship method, slope variation method, and 
variable infiltration capacity (Budyko) hypothesis method. Then, the availability and stability of the three 
methods were compared. The results showed that the runoff in the upper reaches of the Yangtze River 
basin exhibited a decreasing trend from 1956 to 2020, with an abrupt change in 1985. For attribution 
analysis, the runoff series could be divided into two phases, i.e., 1961—1985 (baseline period) and 
1986-2020 (changing period); and it was found that the rainfall-runoff relationship method with 
precipitation as the representative of climate factors had limited usability compared with the other two 
methods, while the slope vatiation and Budyko hypothesis methods had highly consistent results. 
Different factors showed different effects in the sub-basins of the upper reaches of the Yangtze River 
basin. Moreover, human disturbance was the main factor that contributed to the runoff changes, 
accounting for 53.0%—82.0%; and the contribution of climate factors to the runoff change was 
17.0%—47.0%, making it the secondary factor, in which precipitation was the most representative climate 
factor. These results provide insights into how climate and anthropogenic changes synergistically influence 
the runoff of the upper reaches of the Yangtze River basin. 
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1 Introduction 


In recent decades, global warming has changed the hydrological cycle (Hiwasaki et al., 2014; 
Millan, 2014; Poschlod et al., 2020), further impacting runoff via changing precipitation and 
evapotranspiration (Francisco et al., 2022). Furthermore, human activities affect the generation 
and processing of runoff through direct and indirect influences, including reservoir construction 
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and vegetation cover (Zhang et al., 2022; Zhao and Huang, 2022). These two became primary 
driving factors of runoff variability (Ye et al., 2013), clarifying the trends and determining the 
specific contributions of these changes to the runoff amount (Wang et al., 2013; Wang et al., 2022; 
Wu et al, 2022). In the second half of the 20 century, around one-third of 137 globally 
representative rivers with areas ranging from 0.03x106 to 6.30x10° km? experienced about 30.0% 
runoff changes annually, while many of these mid-latitude rivers experienced a 60.0% decrease in 
annual runoff (Walling and Fang, 2003; Blum and Roberts, 2009). These varying degrees of 
runoff variability pose a significant challenge to water resource management. A quantitative 
analysis of the impact of climate change and human disturbance on runoff changes can provide 
theoretical support and decision-making references for the scientific development and utilization 
of water resources (Wang et al., 2013a; Yazdandoost and Morandian, 2021). Most of the existing 
studies investigating the impact of climate change and human activities on runoff used 
hydrological models (Bao et al., 2012), elasticity methods, or empirical statistical methods (Zhai 
and Tao, 2017). These methods can be applied under different conditions. Hydrological models, 
such as soil water assessment tool (SWAT) and variable infiltration capacity (Budyko) models, 
which evaluate the impact of climate changes and human activities by reconstructing the runoff 
time series (Setti et al., 2020), are usually limited by complicated calculation and very large 
amounts of data. Budyko-based elasticity methods use numerous derivative formulas and are 
relatively simple to be calculated (Feng et al., 2017; Liu et al., 2019). Although the physical 
foundations of empirical statistical methods are relatively simple, the satisfying precision and the 
simplicity of calculation process in these methods make them widely used (Hu et al., 2019). 

Previous studies on attribution analyses of runoff changes in arid areas of China were 
extensively studied. Liu and Du (2017) found that human activities were the main factors in 
annual runoff changes in northern and middle-eastern China. Li et al. (2007) indicated that human 
activities accounted for 87.096 of the total reduction in annual runoff in the Wuding River, Loess 
Plateau, China. According to Yuan et al. (2018), climate variability and land use and land cover 
change contributed to 27.1%—49.8% and 50.2%—72.9% decreases of natural runoff in the Yellow 
River between natural and impacted periods, respectively. Du et al. (2011) found that human 
activities and climate changes increased discharge at the rates of 53.0% and 47.0%, respectively. 
The Yangtze River, the world's third-largest river and the main water source for China (Qian et al., 
2012; Peng et al., 2018), provides abundant hydroelectric resources and has become the largest 
reservoir at the global scale. However, the number of natural disasters and water resource 
problems have increased in recent years (Qin et al., 2019; Cao et al., 2021; Hu and Zhang, 2022). 
Consequently, studying the variety of hydrological fluxes under these conditions and investigating 
changes in runoff processes and their impacts on the hydrological cycle of the Yangtze River 
basin attracted the attention of academic community (Li et al., 2016). However, most studies have 
focused on the impact of climate change and human activities on the hydrological processes of the 
entire Yangtze River basin, emphasizing the hydrological changes at the basin outlet. The lack of 
quantitative analysis between different catchment areas makes the comprehensive assessment of 
the impact of climate changes and human activities on the hydrological processes of the upper 
reaches of the Yangtze River basin difficult. Therefore, the focus of this study is on the attribution 
analysis of the changes in water resources in the upper reaches of the Yangtze River basin. The 
changes of runoff in each sub-basin of the upper reaches of the Yangtze River basin were studied, 
and the main causes of these changes were determined. The findings of this study are of the 
utmost importance for the planning and management of water resources in the upper reaches of 
the Yangtze River basin and similar areas around the world. 


2 Materials and methods 


2.1 Study area 
The Yangtze River basin is located in southwestern China (24?27'-35?54'N, 90?33'-112?19'E; Fig. 
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1). The area of the upper reaches of the Yangtze River basin is 100.5x10^ km’, accounting for 
71.4% of the total area of the Yangtze River basin. 
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2.0 Data 


We divided the upper reaches of the Yangtze River basin into several sub-basins according to 
topographical characteristics, including upstream of the Jinsha River (UJR), downstream of the 
Jinsha River (DJR), Mintuo River (MTR), Jialing River (JLR), Wujiang River (WR), and the 
main stream of the upper reaches of the Yangtze River. Eight hydrological stations, including 
Shigu (SG), Pingshan (PS), Gaochang (GC), Fushun (FS), Beibei (BB), Wulong (WL), Cuntan 
(CT), and Yichang (YC) were selected (Fig. 1). The PS station was moved downstream for 2 km 
because of the construction of the Xiangjiaba hydropower station in 2012. The data after 2012 
were collected from the new station, but for the sake of consistency, we keep calling it the PS 
station. In this study, we selected the data from 1956 to 2020. 

Monthly runoff data were obtained from the Hydrological Year Book of China. Meteorological 
data, i.e., mainly monthly precipitation and evaporation, were collected from the China 
Meteorological Data Sharing Service System (http://data.cma.gov.cn). Evaporation data were 
collected from 34 small evaporation dishes in the watershed and calculated using a conversion 
factor of 0.6 between small evaporation dishes and potential real evaporation, which was 
referenced from Sheng (2006). The missing data were interpolated and averaged from 
surrounding areas. Gross domestic product (GDP) and population data in the study area were 
gathered from the Provincial Statistical Yearbooks of Qinghai, Sichuan, Yunnan, and Hubei 
provinces and Chongqing City. Land use datasets, which were developed in the late 1980s and 
2010, were obtained from the Science Data Center of the Chinese Academy of Sciences 
(http://www. resdc.cn). Five images were used in this study, including 1980, 1990, 2000, 2010, 
and 2020. 


2.3 Methods 

2.3.1 Characterization of hydrological variations 

Four methods were used in this study, including the Mann-Kendall (MK) test, cumulative distance 
level method, ordered clustering method, and wavelet analysis. The MK test has good 
applicability in hydrometeorological data analysis because of fewer requirements for sample data 


distribution. The cumulative distance level method accumulates the deviation between measured 
value and mean value (Kong et al., 2016), and visualizes the trend of discrete data points through 
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the curve changes, and then determines the sudden change in precipitation and runoff in different 
periods. The ordered clustering method is used to estimate the transition point of a given time 
series. Finding the critical year is achieved by separating time series into consecutive 
non-overlapping epochs, where the total sum of squared deviations (Sn) before and after is the 
minimum (Zhang et al., 2019). The wavelet analysis is used to identify periodic oscillations of 
signals in standardized runoff, precipitation, evaporation time series, and time-frequency 
correlations. And the Morlet wavelet was chosen for the wavelet transform. 

2.3.0 Attribution analysis 


Climate changes and human activities are generally considered the main factors that impact runoff. 
The rainfall-runoff relationship method, slope changing ratio of cumulative quantity (SCRCQ) 
method, and Budyko hypothesis method were used in this study. 

(1) Rainfall-runoff relationship method 

Because of the correlation between precipitation and runoff, the double mass curve built by 
cumulative runoff and precipitation in the same period can visualize the relationship between 
precipitation and runoff (Pirnia et al., 2019). The main equations are listed as follows: 


Y0-khYP-k,, (1) 
C -50c Os - 4 copy. (2) 
2 AQ O,-Q, 
504 Qabe 199%. (3) 
= AQ Qp-OQ, 
n= 2n. (4) 


AO, 
where Q and P are the average annual runoff (x10* m?) and precipitation (mm), respectively; ki 
and kz are the parameters; AQ is the total change in runoff (x108 m°), Q 4, and Q g are the 


measured runoff in the baseline period and the changing period, respectively (x108 m°); O p ÍS 


the estimated runoff (x108 m°); AQc and AQ y are the amount of climate and anthropogenic 
impacts on runoff, respectively (x10°m?); Cc andC,, are the contribution of climate and human 


activities to runoff change, respectively (%); and ; is the anthropogenic impact benefit. 
Considering different lengths of periods in the base and the changing periods, we performed 
calculations using average values. 

(2) SCRCQ method 

SCRCQ method was essentially used to decompose the effects of climate factors and human 
activities on runoff changes using the cumulative amount, which can avoid the right advertence 
defect from the double mass curve of runoff and precipitation and the multi-point abrupt change 
from the non-parameter MK detection (Wu et al., 2017). Climate impacts in this method include 
precipitation and evapotranspiration. 

The change percentages of runoff (Rsr; 96), precipitation (Rsp; %), and evaporation (Rsz;%) are 
calculated as follows: 


-5S 
Rog = ELS (5) 


Rsp = E x100% , (6) 


WANG Xingbo et al.: Assessment of runoff changes in the sub-basin of the upper... 465 


R= [Sn ioon, (7) 
S ga 

where Sra and Srb are the cumulative runoff in the baseline and changing periods, respectively 
(x108 m?/a); Sp and Sp, are the cumulative precipitation in the baseline and changing periods, 
respectively (mm/a); and Sz, and Sz; are the cumulative evaporation in the baseline and changing 
periods, respectively (mm/a). 

Normally, because precipitation 1s positively correlated with runoff and evaporation is 
negatively correlated with runoff, the contributions of precipitation (Cp; %) and evaporation 
(Cz; %) to runoff are defined as follows: 


C, = SP. 10094, (8) 
SR 
Rsg 
C, =—=x100%. (9) 
SR 
The contribution of climate change (Cc; %) is calculated as follows: 
Coc -CptCg. (10) 


(3) Budyko hypothesis method 
Based on the water balance equation, we described the balance in a basin as follows (Xu et al., 
2014; Wu et al., 2017; He et al., 2019a): 
P=E+R+AS, (11) 
where P is the precipitation (mm); R is the runoff (mm); £ is the actual evapotranspiration (mm); 
and AS is the variation of water storage in the basin (mm), which is considered zero in a closed 
basin over a long period, such as 5-10 a. 
We calculated the long-term annual evapotranspiration based on the following equation 
proposed by Zhang et al. (2001): 


E 
1+o— 


E - d (12) 
1+ PE + e 
P P 

where o is the plant-available water coefficient relating to the way of plants use soil water for 
transpiration (Xue et al., 2022); and Eo is the potential evapotranspiration (mm). Assuming P and 
Eo in Equtaion 12 are independent variables, and we can write Equation 12 as R=f(P, Eo). Taking 
a partial derivative of it, and we can write the total differential of R as follows: 

dR  OR/R dP 4 OR/R dE, 

R P/P P OEy/Ey Ey” 

where d is the total differential; and 9 is the partial derivative. 
The ratio of the growth rates of two interrelated indicators in a given period is defined as the 
coefficient of elasticity, which can easily assess the influence of the variables! change on runoff 
easily (Zhang and He, 2016). The precipitation elasticity of runoff (£p) and the potential 


(13) 


evaporation elasticity of runoff (£x ) are calculated as follows: 


OR/R OR/R 
Ep = RR ; Ey = / E (14) 
9P/P 9E, | Ey 
Then, AQc can be expressed as follows: 
APQ AE,Q 
AQc =€ +E 15 
Qc = £p P Ey E, (15) 


Cc can also be calculated as follows: 
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^ 
Co = ĉc q99s,, (16) 
AQ 


3 Results 


3. Changes in hydrological variables 


The variations in runoff and precipitation in the upper reaches of the Yangtze River basin are 
shown in Figure 2. According to the curve direction, runoff and precipitation exhibited different 
trend changes. The runoff from SG and PS stations fluctuated, the overall trend increased, and the 
change rates were 672.60x10* and 227.00x10* m?/10a, respectively, of which the maximum 
annual runoff from the SG station was 542.00x10* m? in 1965 and the minimum value was 
294.00x 108 m? in 1994. The remaining stations showed a decreasing trend, with the rate shown in 
Table 1. 
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Fig. 2 Annual runoff trends of eight hydrological stations and surface precipitation in six sub-basins of the 
upper reaches of the Yangtze River basin. (a), UJR (upstream of the Jinsha River); (b), DJR (downstream of the 
Jinsha River); (c), MTR (Mintuo River); (d), JLR (Jialing River); (e), WR (Wujiang River); (f), main stream of 
the upper reaches of the Yangtze River basin. SG, Shigu; PS, Pingshan; GC, Gaochang; FS, Fushun; BB, Beibei; 
WL, Wulong; YC, Yichang; CT, Cuntan. The abbreviations are the same in the following figures. 


The MK test was used to determine the significance of trends further, and the values of |Z| 
statistic factor for each station are shown in Table 1. The significance level of the trend test is set 
as 0.05, that is, when |Z|-1.960, the trend of the series change is significant, passing a 95% 
significance test. The maximum test value (1.563) occurred at GC station, and it was less than 
1.960 and failed to pass the significance test. The annual runoff of the remaining seven stations 
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showed a non-significant trend. Among them, the runoff in UJR and DJR sub-basins exhibited a 
non-significant increasing trend, with a rate of 0.675 in the upstream area and 0.223 in the 
downstream area, while the rest of the sub-basins showed a non-significant decreasing trend too. 
For precipitation, the trends of surface precipitation in six sub-basins had different results; only 
WR exhibited a test value greater than 1.960, which passed the 95% significance test. Surface 
precipitation in the other five sub-basins shows a non-significant trend because they did not pass 
the trend test. 


Table 1 Mann-Kendall test results of runoff and precipitation 


: ] 3 Runoff Precipitation 
Station Location Area (km"^) 
Rate IZ Rate IZ 

SG UJR 215,840 0.675 1.551 0.375 1.148 
PS DJR 256,328 0.223 0.724 0.123 0.498 
GC 135,378 —0.935 1.563 

MTR 0.561 0.191 
FS 19,613 —-0.212 1.303 
BB JLR 160,927 —1.495 1.313 —0.443 0.638 
WL WR 87,731 —0.290 0.250 —1.880 2.347 
CT -2.073 0.645 

Main stream 100,837 -1.317 1.887 
YC —1.430 0.408 


Note: SG, Shigu; PS, Pingshan; GC, Gaochang; FS, Fushun; BB, Beibei; WL, Wulong; YC, Yichang; CT, Cuntan; UJR, upstream of the 
Jinsha River; DJR, downstream of the Jinsha River; MTR, Mintuo River; JLR, Jialing River; WR, Wujiang River. The abbreviations are the 
same as in the following tables. |Z], statistic factor value; *, P<0.05 level. 


Cumulative anomaly curves indicating the changing characteristics of runoff at eight stations 
from 1956 to 2020 are shown in Figure 3. More than one inflection point appeared in the 
inter-annual variation of runoff in the upper reaches of the Yangtze River basin. The MK test and 
the ordered clustering method were employed to assess the results further. Taking the BB station 
as an example, images generated by the other two methods are shown in Figure 3b and c. The 
intersection point of two statistic series curves (UF and UB) within the range of the 0.05 
significance level was the mutation point identified by the MK test; similarly, the ordered 
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Fig. 3 Mutation test results of runoff. (a), cumulative anomalies of eight hydrological stations; (b), results of the 
Mann-Kendall test, the red points are intersection points of UF and UB (two statistic series curves); (c), ordered 
clustering method for the BB (Beibei) station, and the red point is the minimum S, (total sum of squared 
deviation). 
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clustering method identified the mutation point based on a very small value. The runoff mutation 
points were examined using three separate methods; diagrams of the two remaining methods for 
the other seven stations are given in detail in Figures S1 and S2, and the results are summarized in 
Table 2. 

It is evident that SG, PS, GC, BB, WL, and CT stations had mutation points near 1985, among 
them, SG, PS, and BB stations showed mutation points identified by multiple methods, suggesting 
that 1985 could be considered as a significant mutation point. By considering factors such as the 
time length and the mutation sites identified by the three methods, and taking the comparability of 
results and process simplification into account, we ultimately identified the year 1985 as the 
mutation point of the entire upstream area. 


Table2 Summary of the three mutation test methods 


Station Cumulative anomaly method Mann-Kendall test Ordered clustering method 
SG 1986*, 1997 1957, 1961, 1967, 1985* 1986*, 1988, 1997, 1960 
PS 1997 1957, 1962, 1968, 1984*, 2012 2005, 1997, 1985* 
GC 1993, 1968 1972, 1986*, 1992, 1993 2017, 1968 
FS 1968 1959, 1966 1966, 2017 
BB 1985* 1960, 1971, 1983* 1985*, 1990, 1993, 1968, 2018 
WL 2002 2004, 2009, 2013, 2019 2019, 2002, 1985*, 1962 
CT 1968, 2017 1957, 1961, 1969, 1986*, 1990, 2018 1968, 2017 
YC 2000, 1968, 1993 1994, 1996, 2000, 2018 2019, 1968 


Note: * represents the mutation point near 1985. 


The Morlet wavelet method was used to calculate the continuous one-dimensional complex 
wavelet transform for the 65-a dataset of measured annual runoff data at each hydrological station. 
The real part of the wavelet transform coefficients provides information about the distribution and 
phase of the signal over time at a specific time scale. Positive values in the graph indicate high 
runoff, corresponding to periods of abundant water, while negative values indicate low runoff, 
indicating dry periods. 

The contour graphs of the real part of the wavelet coefficient of eight stations are shown in 
Figure 4. On the one hand, the runoff time series exhibited strong multi-timescale characteristics 
across all stations. In the first main period, spanning over a cycle of 55—56 a, SG, PS, WL, CT, 
and YC stations went through two transition periods of water abundance and dryness. On the 
other hand, GC, FS, and BB stations had a main period of 43—45 a, representing the most 
significant periodicity observed throughout the entire studied timeframe. For SG and PS 
stations, the second main period is 28 a, characterized by four alternating cycles of positive and 
negative signals within this specific timeframe. GC station exhibited the second main period of 
57 a. FS and BB stations showed the second main period of 12-14 a, with a total of nine 
abundance-depletion transition periods occurring during this cycle. WL station experienced 
approximately three abundance-depletion transition periods within a cycle of 35 a. In CT and 
YC stations, the second main period ranged from 43 to 44 a, encompassing two 
abundance-depletion cycles. However, the second main cycle exhibited relatively weak 
periodicity, localized within specific regions. In annual runoff cycles shorter than 10 a, rapid 
changes occurred more often, characterized by a scattered distribution of extreme value points, 
which suggested that annual runoff at smaller scales fluctuated more frequently and displayed 
evident oscillation behavior. 

Wavelet variance analysis can also be utilized to analyze the periodic intensities. The wavelet 
variance curves and the wavelet coefficient curves of different main periods are illustrated in 
Figure 5. Significant peaks appeared in the wavelet variance of surface precipitation in each of the 
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six sub-basins. The highest peak corresponded to the first main period, representing the periodic 
change in precipitation, while the next smaller peak aligned with the time scale of the second 
main period. Consequently, both runoff and precipitation exhibited similar multi-time scale 
characteristics with strong periodic features at larger time scales, confirming the significant 
influence of precipitation on runoff. 
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Fig. 4 Morlet wavelet coefficients of runoff of eight hydrological stations. (a), PS; (b), SG; (c), GC; (d), FS; (e), 
BB; (£), WL; (g), CT; (h), YC. Red dotted box shows the first main period, and black dotted box shows the second 
main period. 


3.2 Quantitative attribution of runoff variations 
3.2.1 Rainfall-runoff relationship method 


We categorized the precipitation and runoff series of the basin into a base period (1961—1985) and 
a changing period (1986-2020) according to the year of runoff mutation. The double mass curve 
showed the good agreement (Fig. 6). The correlation coefficient of R? values was higher than 0.99, 
and there was no distinct transition point in the double mass curve, which agrees with the linear 
regression curve results. 

Based on the accumulated rainfall-runoff relationship method, climate changes and human 
activities may either negatively or positively influence the total change in runoff. The 
relationships between precipitation and runoff are given in Table 3. Since the total contribution is 
normalized to 1, the sum of the contributions of climate changes and human disturbance must 
exceed 100.095 when they exert opposite effects. Therefore, the signs of the data listed in Table 3 


470 JOURNAL OF ARID LAND 2024 Vol. 16 No. 4 


only indicate the direction of the dominant factor. The results revealed varying contributions of 
precipitation to runoff in different sub-basins. In UJR and main stream, there were opposite 
contributions of climate changes and human activities. The contribution of climate changes was 
111.6% for UJR and 213.9% for main stream. For DJR, JLR, and WR sub-basins, the 
contributions of climate change to runoff changes were 50.896, 61.496, and 70.096, respectively. 
These values indicated that climate changes were the primary factors influencing runoff changes 
in these specific areas. In MTR, the contribution of climate influence was 45.6%. However, the 
main factor affecting runoff changes in this sub-basin was related to human activities. 


= 16 L 1 150 
£ ES 
= 3 -so 
x 127 § 150 
8 = 1960 1970 1980 1990 2000 2010 2020 
S " 5 100 362 
A 1 
3 Ar E 10 
9 1980-1870 1980 1993 2000 3010 30: 139990 20 
. (a) UR (b) DIR 
0 20 40 60 0 20 40 60 
s 60+ E 170 di 
€ t 20 
x oy r- 
o m2 60 0 1980 D 2000 2010 2020 
2 36 L : i € 1980 1990 0 10 202! 
3 zan 110 
E 12 L = eO 1970 ui 1990 2000 2010 202 
S d) WR 
> ; 1 
e 
d. 
3 
8 
S 
> 
2 14 L 1960 1970 1980 1990 2000 2010 2020 
E Year 
$ (e) JLR (f) Main stream 
0 ; i 
° 0 20 40 60 0 20 40 60 


Scale (a) Scale (a) 


Fig. 5 Wavelet variance curves and wavelet coefficient curves of precipitation in different main periods of UJR 
(a), DJR (b), MTR (c), JLR (d), WR (e), and main stream (f) of the upper reaches the Yangtze River basin. Red 
line is the real part curve of the first main period, and blue line corresponds to the second main period. 


3.2.2 SCRCQ method 


The rainfall-runoff relationship method typically focuses on rainfall as a climate factor and does 
not take into account the potential influence of evapotranspiration. However, by introducing the 
SCRCQ method, we considered both rainfall and evapotranspiration simultaneously. In cases 
where runoff variations in a watershed are affected only by changes in rainfall, the slope of the 
cumulative runoff-year relationship should be equal to that of the cumulative rainfall-year 
relationship. Irrespective of other factors influencing the system, the contribution of rainfall to 
runoff changes can be determined by calculating the change ratio of the cumulative rainfall slope 
to the cumulative runoff slope. Similarly, in this study, we divided the study periods into a 
reference period and a changing period using the year 1985 as the dividing point. The calculated 
slopes and change rates for both periods are shown in Table 4. 
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Fig. 6 Double mass curve between cumulative annual precipitation and annual runoff in JLR 


Table 3 Contributions of climate factors and human activities to runoff 


Variation (x 108 m?) Contribution (%) 
, . Measured value Estimated value : 
Sub-basin Period (x10 m?) (X105 m3) Climate Human Climate Human 
factor activity factor activity 
1961-1985 414.70 - = - = = 
UJR 
1986-2020 518.86 530.97 116.27 -12.11 111.60 -11.60 
am 1961-1985 1409.50 z 7 1 E x 
1986-2020 1445.93 1427.90 18.49 17.95 50.80 49.30 
1961-1985 996.20 - - - - - 
MTR 
1986-2020 943.40 972.10 —24.10 -28.70 45.60 54.40 
IR 1961-1985 722.00 - - - - - 
1986-2020 614.84 656.23 —65.77 —41.38 61.40 38.60 
iue 1961-1985 512.20 - - - - - 
1986-2020 471.38 483.62 —28.58 —-12.24 70.00 30.00 
, 1961-1985 4404.60 - - - - - 
Main stream 
1986-2020 4240.43 4053.37 —351.23 187.06 213.90 —113.90 


Note: - denotes a negative effect on the total change. 


Table 4 Slopes in A (runoff), P (precipitation), and E (evaporation), and the results of SCRCQ (slope changing 
ratio of cumulative quantity) 


Sub-basin Rss (96) Rsp(%) Rse(%) C» (96) Cz (94) Cc (96) Cu (96) 
UJR 0.09 0.06 0.02 70.92 —17.08 53.84 46.16 
DJR 0.05 0.02 —0.01 46.76 21.84 68.60 31.40 
MTR —0.05 —0.01 0.01 15.64 21.73 37.37 62.63 
JLR —0.14 —0.04 0.03 26.60 19.76 46.36 53.64 
WR —0.09 —0.03 —0.00 38.20 —].44 36.76 63.24 
Main stream —0.03 —0.01 —0.01 32.78 —15.10 17.68 82.32 


Note: Rsr, Rsp, and Rsg are the change percentages of runoff, precipitation, and evaporation, respectively. Cp and Cz are the 
contributions of precipitation and evaporation to runoff, respectively. Cc and Cy are the contributions of climate changes and human 
activities to runoff, respectively. 


As shown in Table 4, the six sub-basins showed varying slopes in terms of cumulative runoff, 
precipitation, and evaporation due to diverse hydro-meteorological and topographical conditions 
prevailing in the respective sub-basins. Among them, JLR exhibited the highest slope change in 
cumulative runoff of 0.14%. UJR had the highest slope change in cumulative precipitation, with a 
value of 0.06%, and evaporation also showed the highest change at JLR of 0.03%. When 
analyzing the contribution of climate changes to runoff, we found that the sub-basin DJR 
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displayed the highest contribution rate (68.60%), indicating a substantial impact of climate 
changes on runoff, while the main stream showed the lowest contribution rate (17.68%). 

3.2.5 Budyko hypothesis method 

According to the Budyko water-energy balance theory, we divided the study period into the 
baseline period and changing period using the year 1985 as the dividing point. The contribution 
rates of climate variability and human disturbances to runoff changes were quantified, and the 
results are shown in Table 5. 

The mean elasticity coefficients of runoff to precipitation and evapotranspiration are shown in 
Table 5, where e» had a positive value in all catchments and periods. The values ranged from 0.14 
to 0.96, with an average of 0.64 in the baseline period and an average of 0.63 in the changing 
period. Moreover, both the baseline and the changing periods exhibited clear spatial variability in 
the precipitation elasticity coefficient. With an increasing trend in UJR and DJR and decreasing 
trends at the other sub-basins, e» was generally smaller than those of other sub-basins. The ez 
value is negative in all six sub-basins, showing the same spatial trend with ep. Table 5 also shows 
that human activities have a significant impact on the runoff variability in most sub-basins, except 
for DJR. Among them, the main stream exhibited the largest anthropogenic contribution, reaching 
up to 74.3396. 


Table 5 Parameter and result of Budyko hypothesis method 


Aridity index Ep £g 
Sub-basin Baseline Changing Baseline Changing Baseline Changing AQ (m) Cc(%) Cu(%) 
period period period period period period 
UJR 1.99 1.70 0.14 0.20 —0.08 —0.11 8.58 42.08 57.92 
DJR 1.12 1.06 0.47 0.51 —0.26 —0.28 24.43 67.11 32.89 
MTR 0.79 0.83 0.79 0.74 —0.45 —0.42 —22.57 42.74 57.26 
JLR 0.78 0.89 0.80 0.68 —0.46 —0.38 —41.16 38.40 61.60 
WR 0.67 0.68 0.96 0.95 —0.56 —0.55 —14.65 35.91 64.09 
Main stream 0.88 0.89 0.69 0.68 —0.39 —0.38 —69.16 25.67 74.33 


Note: Aridity index can be calculated by Eo/P. ep is the precipitation elasticity, eg is the potential evaporation elasticity, AQ. is the 
amount of climate impacts on runoff; Cc and Cy are the contributions of climate changes and human activities to runoff, respectively. 

In Budyko hypothesis method, the parameter œ (plant-available water coefficient) is 
characterized by uncertainty. The @ value is related to soil and vegetation conditions, primarily 
reflecting the influence of surface characteristics on water balance. It is generally taken within the 
range of 1.0—3.0. In this study, c was set to 2.0 for calculation, and sensitivity was tested within 
the range of 1.0—3.0. The results indicated that the contribution rates of climate changes were 
minimally affected by the selection of œ within this range (Fig. 7). 
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Fig. 7 Parameter sensitivity in Budyko hypothesis method of UJR (I), DJR (II), MTR (I), JLR (IV), WR (V), 
and main stream (VI) of the upper reaches of the Yangtze River basin. Cc is the contribution of climate changes to 
runoff, and o is the plant-available water coefficient. 
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3.2.4 Comparison among the three methods 

Table 6 summarizes the results of the three methods. First, the observed discrepancy can be 
attributed to differences in calculation methods, since the precipitation-runoff relationship method 
normalizes the total contribution rate in the calculation process, summing up to 1. In cases where 
climate and anthropogenic contributions have opposite directions, it employs arithmetic averaging 
to represent the magnitude of changes. This can lead to situations where the absolute contribution 
exceeds 100.0096. Second, the latter two methods consider rainfall and evaporation as 
representative factors of climate changes, demonstrating higher reliability compared with the 
precipitation-runoff relationship method, which only considers precipitation as a factor. Therefore, 
only the results of the other two methods, which exhibit higher comparability and consistency, are 
considered in further comparative analysis. Comparing the results of the other two methods, 
which had high consistency of results, we found that the largest difference in results was only 
11.76% in UJR, while the remaining sub-basins had a mean difference of only 4.60%. 


Table 6 Summary of the results of the three methods 


. Precipitation-runoff relationship SCRCQ Budyko hypothesis method 

iig Co(%) Cu(%) Co(%) Cu(%) Co(%) Cu(%) 
UJR 111.60 —11.60 53.84 46.16 42.08 57.92 
DJR 50.80 49.30 68.12 31.88 67.11 32.89 
MTR 45.60 54.40 37.37 62.63 42.74 57.26 
JLR 61.40 38.60 46.36 53.64 38.40 61.60 
WR 70.00 30.00 36.76 63.24 35.91 64.09 
Main stream 213.90 113.90 17.68 82.32 25.67 74.33 


Note: Cc and Cy are the contributions of climate and human activities to runoff, respectively. SCRCQ, slope changing ratio of 
cumulative quantity. 


3.3 Adaptation and parameter sensitivity analysis 


To assess the adaptability of the attribution methods, we extended the analysis from the mutation 
point in the year 1985 forward and backward to 1982-1988. The results of the attribution analysis 
of the seven change points were analyzed using slope variation method and Budyko hypothesis 
method to compare the relationship between the calculated results of the mutation point in 1985 
and the rest of the points, and determine the adaptation of the two methods. 

Figure 8 shows the contribution rates of climate change to runoff using slope variation method 
and Budyko hypothesis method. The contribution rate of climate change in 1985 fell in the range 
of results for the seven points, while the contribution rates of climate change of the seven points 
in the six sub-basins shown in Figure 8 were consistent. The slope variation method exhibited 
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Fig. 8 Boxplot graph of adaptation of slope variation method (a) and Budyko hypothesis method (b). Cc is the 
contribution of climate factors to runoff. Boxes indicate the IQR (interquartile range, 75" to 25" of the data). The 
median value is shown as a line within the box; square is shown as mean; and black diamond is outlier. Whiskers 
extend to the most extreme value within 1.5xIQR. 
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better adaptation than Budyko hypothesis method. Within the same basin among seven points, 
DJR exhibited the largest range of variation, still not exceeding 20%, while JLR and main stream 
varied only 396—496. However, compared with slope variation method, Budyko hypothesis 
method showed an anomaly in JLR, where the climate change contribution value in one of the 
seven years exceeded the normal range of variation. The basin's Cc normally fluctuated between 
32% and 38% during 1983-1990, while the contribution rate in 1982 was lower (23%). These 
findings were related to the larger AZo due to the change in potential evapotranspiration. The Cc 
value in the remaining sub-basins varied mostly in WR over the 7 a, reaching 17%, and UJR had 
the smallest variation range (<6%). 


4 Discussion 


In the second half of the 20 century, around one-third of among 137 globally representative 
rivers with basin areas ranging from 0.03108 to 6.30x10* km? experienced annual runoff changes 
exceeding 30% (Milliman et al., 2008; He et al., 2019b), and many of these mid-latitude rivers 
underwent a decrease in annual runoff by even 60%. While, as for the dominant factor on runoff 
change, human activities including dam construction, land use changes, social development, 
irrigation, and inter-basin water transfers, have been confirmed to be the primary factors in the 
Yellow River and its tributaries (Wang et al., 2012), the Zhujiang River (Wu et al., 2012), and the 
Huaihe River (Liu and Xia, 2013). 


4.1 Impact of climate changes 


Previous studies have shown that global warming-related intensification of the hydrological cycle 
affects runoff by increasing evapotranspiration (Chien et al., 2013; Wang et al., 2016). We found 
that precipitation and runoff in the study area exhibited parallel trends and cyclical changes, 
indicating a strong positive correlation. Taking JLR as an example, we analyzed the runoff from 
BB station and the surface precipitation in JLR (Fig. 9). Precipitation and runoff regression curve 
showed a good linear relationship, and R? of baseline period was smaller than that of changing 
period, and the correlation between precipitation and runoff was reduced by human activities. 
Further details regarding the correlation at the other five sub-basins are given in Figure S3. 
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Fig. 9 Relationship between precipitation and runoff in different periods. 1961—1985, baseline period; 
1886-2020, changing period. 


4.2 Impact of hydraulic engineering 


In addition, our studies implied that large-scale water projects played an important role in runoff 
redistribution, as confirmed in previous study (Zhang et al., 2021). Chinese government 
implemented a series of comprehensive water conservancy hubs with good regulation capacity 
and large reservoir capacity based on the joint operation of reservoir groups in the upper reaches 
of the Yangtze River basin (Xu et al., 2017), including the Jinsha River middle reaches reservoir 
group consisting of six reservoirs, the Jinsha River downstream reservoir group, the Yalong River 
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reservoir group, the Minhe River reservoir group, the Jialing River reservoir group, the Wujiang 
River reservoir group, and the mainstream Three Gorges-Gezhouba cascade reservoir group. The 
total regulation capacity reaches 110.0x10? m?. Furthermore, about 34.0x10? km of dikes and 
more than 40 flood storage areas with a capacity of around 62.7x10? m? also interact with 
reservoirs to influence runoff. About 96 large-scale reservoirs existed by the end of 2020, 
including the Three Gorges-Gezhouba in the upper reaches of the Yangtze River basin, and the 
Liyuan, Xiluodu, and Xiangjiaba dams on the Jinsha River, with a total water storage capacity of 
128.0x10? m?. The reservoirs in JSR have a storage capacity of 648.0x10° m°, the Three Gorges 
reservoir has a storage capacity of 393.0x10* m’, and the group of reservoirs in WR has a storage 
capacity of 111.0x10* m?. Compared with the year 2000, the number of large reservoirs increased 
by 78.096, and the significant increase in the proportion of human activities during this period was 
closely related to the substantial completion of water projects. 


4.3 Impact of land use condition 


The land use condition affecting the hydrological and sedimentary changes in the basin by 
influencing vegetation greening and subsurface texture (Zeng et al., 2016; Mankin et al., 2019). 
For example, the reduction of farmland and the expansion of urban land decreasing soil 
precipitation infiltration rate and increasing runoff (Wang et al., 2023b). The expansion of water 
body can decrease runoff by enhancing evaporation. Also, the larger canopy area of forest land 
strongly enhances the ability of precipitation interception (Raz-Yaseef et al., 2010; Yu et al., 
2022). Based on remote sensing monitoring data on the land use status, researchers constructed 
the land use transfer diagram. The increase in construction land was mainly transferred from 
farmland and forest land, with a significant rise of 231.6%. Also, a 15.8% decrease in unused land 
appeared, mainly shifting toward water body and grassland. These data are consistent with Wang 
et al. (2013b), who found that the Yangtze River basin has experienced a rapid urbanization 
process characterized by farmland reduction and urban land expansion. 


4.4 Impact of socio-economic factors 


Since the mid-20" century, the Yangtze River, as one of the Chinese crucial economic belts, 
experienced robust economic development and a rapid increase in population. The gross domestic 
product and population in the main provinces of the upper reaches of the Yangtze River exhibited 
strong upward trends, which had significantly increased water consumption for industrial, 
agricultural, and domestic purposes, further reducing the area of runoff (Wang et al., 2013b). 
Moreover, a land reform happened around 1980, farmers were encouraged to manage their lands, 
and water demand for agriculture irrigation was increased (Yang and Tian, 2009). According to 
the Bulletin of Water Resources in China, water consumption for the Yangtze River was 
1737.1x10? m? in 1997 and increased up to 1957.6x105 m? in 2020. 


4.5 Summary of change in sub-basins 


The Jinsha River exhibits significant variations in topography, background conditions, and other 
characteristics that differ between UJR and DJR (Feng et al., 2017; Liu et al., 2017). Taken as a 
whole, climate changes and human activities influences varied around 50.0% throughout the 
Jinsha River. Also, Zhang et al. (2019) considered that climate changes explained 56.3% of the 
reduction in runoff, which the same conclusion that the runoff influence from climate and human 
activities were considered to be evenly matched (each was around 50.096). MTR is the most 
important water system in Sichuan Province, and the Minjiang River basin covers multiple cities, 
representing a population and industrial agglomeration area (Zhou et al., 2008; Tang et al., 2013). 
However, a significant portion of the basin lies in a water-deficient region and the water shortage 
is addressed by water diversion from the Minjiang River. Such human activities affect runoff to a 
great extent, being the primary influencing factor for the runoff change, accounting for 60.0% 
based on our findings, which is consistent with Guo et al. (2022), being a 69.0% proportion based 
on Budyko hypothesis method. According to our results, the main influencing factor of runoff 
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variability in JLR is human activities, with a contribution rate of around 50.096—60.096, which is 
due to the variety of areas interlaced with mountainous terrain and basins (Zeng et al., 2015; Shao 
et al., 2021), and the large scale of the hydraulic engineering with a total storage capacity of 
78.0x10* m?. WR sub-basin is located in a typical karst mountainous region with a dense 
population and abundant land resources (Wang et al., 2023a), with 11 levels of multi-level 
hydropower schemes ready to be implemented (Fu et al., 2023). There are more than 120 
constructed or constructing hydroelectric stations, with capacities exceeding 500 kW. These 
extensive human activities contribute to approximately 63.0% of the runoff variability in the basin. 
The mainstream of the upper reaches of the Yangtze River basin is the most economically 
developed and densely populated region (Zhang et al., 2006; Xiao et al., 2018), with the primary 
human activities contributing to over 70.0% of the runoff variability, making it the most 
significant factor influencing the basin's runoff changes. Our results are consistent with the runoff 
change attribution results obtained in the study of Zhai and Tao (2017). 


4.6 Uncertainty analysis 


The amount of data is relatively limited due to the data from hydrological stations at the outlet of 
each sub-basin, which may not be sufficiently precise to represent the situation in the entire river 
basin. Additionally, our study treated climate changes and human activities as two independent 
variables for attribution analysis, while they cannot be completely separated and may influence 
each other. Interactions between human activities and climate changes are complex and 
challenging to determine, so further research is needed for their accurate quantification. 


5 Conclusions 


In this paper, we characterized the long-term runoff sequences at major stations in the upper 
reaches of the Yangtze River basin, China. The runoff showed a non-significant decreasing trend 
in other stations except for a slight increase in PS and SG stations of the Jinsha River. The abrupt 
change occurred in 1985. When comparing the results of runoff attribution using the 
precipitation-runoff relationship method, slope change method, and Budyko hypothesis method, it 
was found that the latter two methods exhibited a high degree of consistency and passed the 
adaptability test. Except that JSR showed a different result with a proportion of basically 50.0% 
each for climatic factors and human activities. Human activities in MTR, JLR, WR, and main 
stream sub-basins accounted for from 53.0% to 82.0% of the impact on runoff and were the main 
factors driving runoff change. The contribution of climate changes ranged from 17.0% to 46.0%, 
making them the secondary factor. Finally, climate and various human activities have diverse 
effects under different conditions. Managing according to the characteristics of each sub-basin 
ensures the sustainable use of water resources in the Yangtze River basin, facilitating better 
realization of complementary feedback among the water resources of sub-basins and providing a 
foundation for exploring differentiated water management strategies at the sub-regional scale. 
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Fig. S1 Mann-Kendall mutation test results of runoff for eight hydrological stations. (a), SG (Shigu); (b), PS 


(Pingshan); (c), GC (Gaochang); (d), FS (Fushun); (e), BB (Beibei); (f), WL (Wulong); (g), CT (Cuntan); (h), YC 
(Yichang). The abbreviations are the same in the following figure. 
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Fig. S2 Runoff ordered clustering test results for eight hydrological stations. (a), SG; (b), PS;(c), GC; (d), FS; 
(e), BB; (©, WL; (g), CT; (h), YC. S, is the total sum of squared deviation. 
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Fig. S3 Relationship between precipitation and runoff in different periods for five sub-basins. (a), UJR 
(upstream of the Jinsha River); (b), DJR (downstream of the Jinsha River); (c), MTR (Mintuo River); (d), WR, 
Wujiang River; (e), main stream of the upper reaches of the Yangtze River basin. 1961—1985, baseline period; 


1886-2020, changing period. 


